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Abstract 

Cosmological consequences of adding the Cold Dark Matter (CDM) to the exactly 
solvable stringy Dark Energy (DE) model are investigated. The model is motivated 
by the consideration of our Universe as a slowly decaying D3-brane. The decay of this 
D-brane is described in the String Field Theory framework. Stability conditions of 
the exact solution with respect to small fluctuations of the initial value of the CDM 
energy density are found. Solutions with large initial value of the CDM energy density 
attracted by the exact solution without CDM are constructed numerically. In contrast 
to the ACDM model the Hubble parameter in the model is not a monotonic function 
of time. For specific initial data the DE state parameter woe is also not monotonic 
function of time. For these cases there are two separate regions of time where woe 
being less than —1 is close to —1. 
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1 Introduction 



Nowadays strings and D-branes found cosmological applications related with the cosmolog- 
ical acceleration pQ-j3]- The combined analysis of the type la supernovae, galaxy clusters 
measurements and WMAP data provides an evidence for the accelerated cosmic expan- 
sion jUElini- The cosmological acceleration strongly indicates that the present day Universe 
is dominated by smoothly distributed slowly varying Dark Energy (DE) component (see [7] 
for reviews), for which the state parameter wde is negative 1 . Contemporary experiments 
give strong support that currently the state parameter wde is close to —1, wde = — 1 ±0.1 

ei ei maim. 

From the theoretical point of view the specified domain of w covers three essentially 
different cases: w > — 1, w = — 1 and w < — 1 (see [T2], and references therein). The 
most exciting possibility would be the case w < — 1 corresponding to the so called phantom 
dominated Universe. In phenomenological models describing this case the weak energy 
conditions g>0,g + p>0 are violated and there are problems with stability at classical 
and quantum levels [13]. Thus, a phantom becomes a great challenge for the theory while 
its presence according to the supernovae data is not excluded. 

A possible way to evade the stability problem for a phantom model is to yield the phantom 
as an effective model of a more fundamental theory which has no such problems at all. It has 
been shown in [3] that such a model does appear in the string theory framework. This DE 
model assumes that our Universe is a slowly decaying D3-brane which dynamics is described 
by the tachyonic mode of the string field theory (SFT). The notable feature of the SFT 
description of the tachyon dynamics is a non-local polynomial interaction J1]"EB]- It turns 
out the string tachyon behavior is effectively described by a scalar field with a negative 
kinetic term (phantom) however due to the string theory origin the model is stable at large 
times. 

In ^2] we have found an exactly solvable Stringy DE model in the Friedmann Universe. 
This model is a modified version of the effective SFT model [3] and is inspired by SuperSFT 
calculations J7] . First level calculations in the SFT give fourth order polynomial interaction. 
Higher levels increase a power of the interaction. Exactly solvable model has a particular 
six order polynomial interaction potential. However, small fluctuations of coefficients in 
that potential do not change the solution qualitatively and one can say that the model J2] 
represents the behavior of nonBPS D3 brane in the the Friedmann Universe rather well. It 
is interesting to investigate the dynamics of the model in the presence of the Dark Matter. 
This is a subject of the present paper. 

It turns out from the observational data that DE forms about 73% and the Dark Mat- 
ter forms about 23% of our Universe. Thus because of a significance of the Dark Matter 
component in the Universe in the present paper we investigate an interaction of the phan- 
tom matter considered in ^2] with the CDM. It seems impossible to find exact solutions in 
the presence of the CDM, except the case when the DE state parameter is a constant |19j . 
so we use numeric methods to analyze the behavior of the phantom field and cosmological 
parameters in our model. 

^ere wde is usual notation for the pressure to energy ratio. 
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2 Exactly solvable Phantom Model 

We start by recalling the main facts related to the model considered in [T2] . This is a 
model of Einstein gravity interacting with a single phantom scalar field in the spatially flat 
Friedmann Universe. Since the phantom field comes from the string field theory the string 
mass M s and a dimensionless open string coupling constant g Q emerges. The action is 



>=;(■ 



Ml „ 1/1 



s = / d 4 xV=g + -r +^d^d v <p - v((j>) , (l 



\2Mf~ ' 9 2 V 2 " v 
where Mp is the reduced Planck mass, g^ v is a spatially flat Friedmann metric 

ds 2 = — dt 2 + a 2 (t)(dx\ + dx\ + cfx 2 .). 

and coordinates (t, Xi) and field are dimensionless. Hereafter we use the dimensionless 
parameter m p for short: 

m 2 = (2) 

s 

If the scalar field depends only on time, i.e. <fi = <f>(t), then independent equations of motion 
are 

3# 2 = Qde, Qde = -^0 2 + V((f>), 
mi 2 

(3) 

3H 2 + 2H = - i- p DS = -^U 2 - 7(0), 

Here dot denotes the time derivative, H = a(t)/a(t), Qde and pde are energy and pressure 
densities of the DE respectively. One can recast the system (J3J) to the following form 

i 

H 



2ml 

(4) 

iH " - k {-\^ + v ^ 

Besides of this there is an equation of motion for the field <ft which is in fact a consequence 
of system (jSJ). 

Following the superpotential method [201 (see also [21]) we assume that H(t) is a function 
(named as superpotential) of (f>(t): 

H(t) = W{(j>{t)). (5) 

This still does not give a systematic way to find general solutions to the system (J3J) but 
allows one to construct W((j)) and V((f)) for a known function <fi(t). We take for 4>(t) 

(j>(t) = Ataah(ut). (6) 
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This function is known to describe effectively the late time behavior of the tachyon in the 
4- dimensional flat case [221 I2H] • The function <p(t) satisfies the following equation 



Hence, we obtain 



and corresponding potential 



V 



^ ) = £(*-^' + I ^ (3 ,8) 

We have omitted an integration constant in (|7J) to yield an even potential (jSJ). It is typical 
that to keep the form of solutions to the scalar field equation in the presence of Friedmann 
metric one has to modify the potential adding a term proportional to the inverse of the 
reduced Planck mass Mj, [T21I24] . 

The described solution leads to a number of cosmological consequences. The Hubble 
parameter 

H = ^Ar tanhfwt) (1 - -(tanh(iut)) 2 ) (9) 
2m 2 V 3 J 

goes asymptotically to uA 2 / (3m 2 ) when t goes to infinity. Once H(t) is known one readily 
obtains the scale factor 

/ w »^ /M 2 (cosh(cjt) 2 - 1)\ . . 

.(«) = ao(coshM))^ exp ( ^JJ^' ), (10) 

where ao is an arbitrary constant, and the deceleration parameter 

aa 18m 2 (cosh(wt)) 2 

gm = ~ = — 1 ? . (11) 

« 2 A 2 ({cosh{u t) f - 1) (2 (cosh(c^)) 2 + 1) 

It follows from formula (fTT|) that the Universe in this scenario is accelerating. 
The expression for the state parameter is the following 



^44 =-1-12^4^7 02)2 

Qde(</>) p 2 (3A 2 



w DE {4>) = —r^ = -1 - 12m (12) 



Point <p = A corresponds to an infinite future and therefore wde — 1 ast^oo. 

Plots for the Hubble, deceleration and state parameters are drawn in Fig^ (Hereafter we 
assume A = u = 1 for all plots). 

Thus, we conclude by noting that in our model the phantom field provides the DE 
dominated accelerating Universe. 
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Figure 1: The time evoution of the Hubble parameter H(t) (left), the deceleration parameter 
q(t) (middle) and the state parameter Wde{£) (right) in the exactly solvable model for m 2 v = 
0.2. 

3 Interaction with Cold Dark Matter 
3.1 The Model 

Now we are going to couple in a minimal way a pressureless matter of energy density qm 
(the CDM) to our model such that the Friedmann equations get an extra term 



3H 2 
H 



i 



1 ■ 



<fi 2 - Q M {a)) . 



and the equation describing the evolution of scalar field has previous form 

4> + 3H<p - = 0. 

From (|T3|) - (|T5|) we obtain the conservation of the energy density for the CDM: 

q m + 3Hg M = 0, 

that after integration gives 



L>.M < !>M.»[- 

a 



(13) 
(14) 

(15) 
(16) 

(17) 



where constants Qm,o and a are initial values of Qm and a correspondingly. From (fTTj) we 
obtain equation (|T3*j) in the following form: 



3# 2 = ~ 2 [v{4>) - 4- + QM fl 

mi \ 2 V a 



Following the lines of jTzj we address to our analysis the questions of cosmological evo- 
lution and stability. 
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The straightforward way to study a stability of solutions to the system of equations 
(jT3*j) - (fT5j ) is to exclude qm from (JT3J) , (U3J) and obtain the following system: 



+ 3H<j)-V' = 0, (19) 



1 ' l l2 



2H + 3H* = —^ + V((j>)). (20) 

Depending on the initial values of H, <p and <fi, which have been considered as independent, 
this system describes our model either with or without the CDM. In particular, the initial 
values: H = 0, <p Q = and <p = Au correspond to the exact solution JBJ). 

Either one can perform calculations using nonautonomous system of equations [22] which 
can be obtained from equations (|15|) and (J 18)) in the following form 

# 1 (21) 



dn H((p,ip,n 

# = _3^ + _ 
dn H(<p,ip,n 



dip 1 dV(d>) , s 

3 ^ + WU i \ — IT 1 ! ( 22 ) 



where 



H((p,ip,n) = \j~\^ 2 + + 0M,oe~ 3n , (f) = ip and n = ln(o/a ) 



3.2 Stability analysis for small fluctuations. 

In ^2] we have analyzed the system (|13| ) - (fT5j) without the CDM under condition A = uo = 1 
and found that the exact solution <ft = tanh(t) is stable with respect to small fluctuations of 
the initial conditions if and only if m 2 , ^1/2. 

Let us consider the behavior of the solution of system (fT9"j) - (j20j) in the neighborhood of 
the exact solution 

0o(£) = tanh(t), 

<j> (t) = Mt) = 1 - tanh(t) 2 , 

H (t) = ^2 tanh(t) (l - itanh(i) 2 ^ . 

Substituting 

H(t) = H (t) + eH x (t), (j)(t) = o (t) + e0i(t) and 0(t) = ip(t) = ip (t) + e^t), (23) 



6 



in (|19|) and (|2(J|) we obtain in the first order of e the following equations: 



0i = if>i, 

If „ 6m 2 - 1 \ 



(2 cosh(t) 2 + 1) tanh(t) ^ 
2m 2 cosh(t) 2 1 cosh(t) 2 1 ' 

*• ^ i3F( 1 + ,! - ta ' ) ™ h,1), " 1+ 

1 (1 + 2cosh(t) 2 )tanh(t) 

+ 2cosh(t) 2 ^ 1 2cosh(t) 2 1- 

System (PH) has been solved with the help of the computer algebra system Maple. The 
exact dependence <j>i(t), ipi(t) and H{t) are too cumbersome to be presented here. The main 
result is that for m 2 ^1/2 functions (f>i(t), 4>i(t) and H\(t) are bounded functions and our 
exact solution is stable. 

Note that numerical calculations show that if ml ^1/2 then even for large initial values 
of the CDM energy density numerical solutions tend to the exact solution as t tends to 
infinity. 



3.3 Numeric solutions. Time dependence. 

At this point we pass to numeric methods because it seems impossible to find exact solutions 
in the presence of the CDM. To analyze the cosmological evolution it is instructive to plot 
phase curves for the scalar field as well as evolution of the state parameter wde for the scalar 
matter. In addition we find numerically a ratio of the energy densities for the CDM and the 
DE. Experimental bounds for this ratio is known and estimated to be near 1/3 so we can 
find the time point we live and a corresponding value of wde in our approach. 

Due to equation (|18|) initial data <po, 0o an d H do fix an initial value of the CDM 
density. To have a given initial energy density of the CDM we take 0o an d <fio and find 
the corresponding value H . In particular, to have qm,q — 1, 4>o = 0, O = 1 for m 2 = 0.2 
we must take Hq = a/5/3 ~ 1.29. For this initial values numeric solutions are presented 
graphically in FigJ21 In FigJSl we present the same plots for qm,q = 100. Corresponding 
H = 10a/5/3. Comparing FigsO and 01 with Fig^ one can see that our solutions with 
and without the CDM are different only in the beginning of the evolution where the CDM 
dominates (if exists). Note, that the behavior of the Hubble parameter in the presence of 
the CDM is not monotonic and the DE state parameter may be not monotonic as well. 



3.4 Numeric solutions, (^-dependence. 

It turns out that values of the Wde as well as ratio Qcdm/qde which are observational 
cosmological parameters 00111] can be found easier using equations (jUJ), (}2"2"j) as functions 
of the e-folding number n. However, it is more instructive to find a dependence on and 
not on n. 
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Figure 2: The time evolution of the scalar field 4>(t), the Hubble parameter H(t) and the 
state parameter Wde{^)- 0m,o — 1> m f> — 0.2 and O = 1. 




1 2 3 4 5^6 1 2 3 4 5^6 



Figure 3: The time evolution of the scalar field 0(£), the Hubble parameter H(t) and the 
state parameter WDE(t). qm,o = 100, = 0.2 and 0o = 1- 

Let us recall that from an analysis of our phantom model without the CDM we know 
that the scalar field interpolates between an unstable and a nonperturbative vacua during 
infinite time similar to the non-BPS string tachyon [221 EHj- In our notations nonperturbative 
vacuum corresponds to = +1. In the pure phantom model the evolution is described by 
4>(t) = Atanh(ujt) function, where A and to can be rescaled to 1. This dependence is 
monotonic and this allows us to find physical variables as functions of 0. 

The situation is more complicated in the presence of the CDM. First, we do not know 
an exact time dependence of the scalar field. Second, it is not evident for arbitrary initial 
data and value of parameter w? v that the scalar field evolves monotonically. However, in 
the particular cases presented in Figs. 121 and El our solutions 0(t) are monotonic functions 
of time and moreover look like tanh(t) at large times. Below we are interesting in solutions 
which approach the nonperturbative vacuum during an infinite time. Thus, the point = 1 
corresponds to an infinite future. The 0(t) dependence can be found numerically to pass 
from coordinate to the time. 

In FigsEHHlwe plot results of numeric solutions to equations (j2U), fj2^|) that allow us to 
find physical variables such as H, wde, Qcdm/qde as function of the field 0. These sets of 
plots differ in an initial velocity of the scalar field. Note that it follows from (|18j) that there 
exists a maximal initial velocity tpom for our phantom field. tpo m depends on values of 0o, 
qm,o and ao and does not depend on m^. In all plots a = 1 and 0o = . All plots have 
three curves: black ones are phase curves, red ones are w-s and blue ones are Qcdm/qde 
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= 0.01, 



Pm,o = 1 



Pm,o = 100 




Figure 4: ^-dependence of the velocity </> (black line), the state parameter wde (red line) 
and the Qcdm/qde ratio (blue line). Initial velocity of the scalar field is equal to 1 and 
m 2 p = 0.2. 

ratios. In Fig. 01 the initial velocity is equal to 1 (which is the same as for the exact solution), 
ml = 0.2 and qm,o is equal to 0.01, 1 and 100 from left to right. Here we see that the scalar 
field reaches +1. This indicates a stability of the system with respect to fluctuations of the 
initial CDM energy density for small m 2 v . In FigJSJthe initial velocity is equal to 0.72. The 
first row there corresponds to m 2 p = 0.2 and qm,o is equal to 0.01, 1 and 100 from left to 
right. The second row shows the behavior of the system with equal to 0.6 and 1 and 
Qm,o equal to 0.01 and with equal to 1 and qm,o equal to 1. One again sees from these 
plots that the scalar field reaches 1 for small values of ml in a wide range of an initial CDM 
energy density. 

This situation is broken for greater even for a small initial CDM energy density and 
the field does not reach 1. In Fig. E1£?m,o is taken to be 1, the initial velocity is equal to 
its maximal possible value ipom and ml is equal to 0.2, 0.6 and 1 from left to right. Here 
we again observe a stability for small m p and also find out that for large the scalar field 
goes beyond the point 1. Also for the maximal possible initial velocities wde and qcdm/ Qde 
functions have a discontinuity. One can understand this qualitatively because the energy 
density for the DE has two terms with opposite signs. Indeed, the scalar field is a phantom 
and its kinetic energy is negative while the potential term is positive. Thus at some point 
the energy density of the DE changes the sign and develops a discontinuity in w^e and 
Qcdm I Qde ratio. Such a behavior is rather undesirable from cosmological point of view, 
since the is no observational data indicated singular behavior of cosmological parameters 
and we do not consider corresponding plots further. 

Hence, seeking for a situation where field <p approaches 1 and there is no cosmological 
singularities during this evolution we are left with the first row in Fig. 0] and Fig. E3 In 
this plots phase curves show that field indeed depends monotonically on time because <fi 
is always positive during the evolution. Looking for specified plots we draw the reader's 
attention to the following interesting properties of our model. First, Qcdm/ Qde ratio depen- 
dence is monotonic and experimentally measured value 1/3 is close to the beginning of the 
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m 2 p = 0.2, p M)0 = 0.01 m 2 p = 0.2, p Mj0 = 1 = 0.2, p M , = 100 




ml = 0.6, pa/,0 = 0.01 m 2 v = 1, p Mj0 = 0.01 mjj = 1, p Mfi = 1 




Figure 5: 0-dependence of the velocity (ft (black line), the state parameter wde (red line) 
and the Qcdm/qde ratio (blue line). Initial velocity of the scalar field is equal to 0.72. 

evolution. For example, in Fig. |^ (left) this point corresponds to wde ~ —1-02 and <fi w 0.09. 
Moreover, this is not a distinguished value and it follows from the model that this ratio will 
decrease with time. Second, for pm,o large enough wde behaves non-monotonically. 

4 Discussion and Conclusion 

To summarize, let us also note that we get an existence of a region of the initial energy 
density of the CDM, for which Wde is not monotonic. Such a behavior is interesting and 
very surprising. We see that for large initial energy densities of the CDM wde grows with 
time from minus infinity to approximately —1, then goes down to a local minimum and after 
this grows again asymptotically approaching —1. Note, that it has been proved in !'2(J that 
under some (compare with j27j) conditions the phantom matter cannot cross the w — — 1 
barrier. For small initial energy densities of the CDM we cannot say that wde has a local 
minimum but its rate of change becomes slower and faster. We especially point out that the 
very first time stamp when wde approaches —1 (or its rate of change decreases essentially) 
is approximately the same where Qcdm/ ' Qde ratio is close to 1/3. It is worth to note that 
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Figure 6: 0-dependence of the velocity (black line), the state parameter wde (red line) 
and the Qcdm/qde ratio (blue line). Initial velocity of the scalar field is equal to its maximal 
possible value. 



juxtaposing Figures I2H21 and Figures EHH1 we see that the present ratio of the CDM and the 
DE energy densities also corresponds to a local minimum of the Hubble parameter. 

It is interesting to find an influence of the higher open string mass levels as well as an 
influence of the closed string excitations on the obtained picture. Even in the flat space-time 
the dynamics of a D-brane change drastically when the closed string excitations are included 

mi 121. 

We get a stable behavior and smooth cosmological parameters in the stringy inspired 
model only in the case when the dimensionless parameter m 2 is less than 0.5. This restricts 
the parameters of the original theory [3] the model considered in this paper comes from. Let 
us recall that m 2 , is related with the reduced Planck mass, the string mass parameter M 2 
and the open string coupling constant g%: 

n 2 M 2 
,2 _ iJo 1V1 P 



Therefore, to have an acceptable cosmological solution we have to assume that 2g^Mp < M 2 . 
Here Mp is the 4-dimensional Planck mass. The effective string mass depends on the 10- 
dimensional mass parameter l/l s , where l s is a string length and the compactification volume 
Vq. Taking into account that M 2 ~ Vq/1 8 and using the usual requirement that l s ~ 1/M P 
we get 

2 = 9 M P g M P L s gj^ 
p M 2 V 6 V 6 <U ^' 

that looks reasonable from a general point of view. 
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